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Abstract 

The subject of this study is an iterative Bermudan option pricing algorithm based on 
(high-dimensional) cubature. We show that the sequence of Bermudan prices (as 
functions of the underlying assets' logarithmic start prices) resulting from the itera¬ 
tion is bounded and increases monotonely to the approximate perpetual Bermudan 
option price; the convergence is linear in the supremum norm with the discount fac¬ 
tor being the convergence factor. Furthermore, we prove a characterisation of this 
approximated perpetual Bermudan price as the smallest fixed point of the iteration 


procedure. 



When Nicolas Victoir studied “asymmetric cubature formulae with few points” j2] 
for symmetric measures such as the Gaussian measure, the idea of (non-perpetual) 
Bermudan option pricing via cubature in the log-price space was born. In the 
following, we will discuss the soundness and convergence rate of this approach 
when used to price perpetual Bermudan options. 

Consider a convex combination (oi,..., am) S [0, (that is, J2T=i = 1) 
and xi,...,Xm £ Then there is a canonical weighted arithmetic average 
operator A associated with a, x given by 

m 

k=l 

Now suppose c G (0,1), g,h Ag > g, Ah = h and OW g < h. Define an 

operator V on the cone of nonnegative measurable functions by 

V : f ^ {c- Af)V g. 

Since A is positive and linear, thus monotone (in the sense that for all /o < /i, 
^/o < ^/i). it follows that V must be monotone as well. Furthermore, whenever 
> /. we have that AVf > 'Df, as the linearity and positivity of A combined 
with our assumption on g imply 

AVf > cA^f W g>{c- Af) V g = Vf. 

Finally, due to our assumptions on h and g, we have for all nonnegative f < h, 

Vf < cAh < h. 

Summarising this, we are entitled to state 

Lemma 1. Adopting the previous paragraph’s notation and setting 

Q:={f<h : Af>h}, 
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we have that 


V-.Q^Q, 

V is monotone (i e order-preserving), and AD — V is nonnegative. 

This is sufficient to prove 
Theorem 1. For all n S Nq, 

D^+\gV0)>V^igW0)=:qn. (1) 

Furthermore, 

q ■.= lim V'^{g V 0) = sup T»”(g V 0) S Q 
n6No 

and q is the smallest nonnegative fixed point ofD. 

Proof. 1. The proof of equation o is a straightforward induction on n where 
we have to use the monotonicity of D in the induction step. 

2. Since D maps Q itself, the whole sequence {D'^{g V is bounded by 

h. This entails q < h as well. Using the linearity of sup and our previous 
observation that AD — D>0 (Lemma we can show 

VyeR'" Aq{y) = sup ^ (P"(g V 0)) (y) 

n^No 

> sup T >”(5 V 0)(y) = q(?/), 

n^No 

which means Aq > g. As we have already seen, q < h, so q G Q. 

3. Again, due to the linearity of sup, D and sup„gpj^ commute for bounded 
monotonely increasing sequences of functions. Thereby 

Dq = sup DD^{g V 0) = supT >"(5 V 0) = g. 

nGNo nGN 
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4. Any nonnegative fixed point p of P must be greater or equal gWO. There¬ 


fore by the monotonicity of sup and V, 

sup P"p > sup P"(p V 0) = q. 
rtSNo nSNo 


□ 


Lemma 2. Using the previous Theorem’s notation, we have for all x and 
nSNo, if qn+i{x) = g(x), then qn{x) = g{x). 

Proof. By the monotonicity of the sequence (gn)nGNo (Theorem^, we have 
g{x) < qoix) < qnix) < qn+lix). 


□ 


Theorem 2. For all n S N, 

||<?ra+l ~ ^ ' \\Qn ~ qn-1 llc‘>(R‘^,R) ' 

Proof. The preceding Lemma [21 yields 

||<Zn-|-l — <ZTt|lcO(R<i R) = Ikn+l “ 

= l|c ■ Aq„ - ((c ■ Aqn-i) V 5)IIco({c.A9„>3}.r) 

via the definition of qt+i as (cAqi) V g for i = n and i = n -|- 1. But the last 
equality implies 

Ikn-l-l — 9n|lcO(R<i,R) ^ II’ ^9™ — C • A(7„_i || >g},R) 

< ||c • Aqn — C ■ Aqn-l llcOj-Rd^R) ■ 

Since A is linear as well as an L°“-contraction (and therefore a C^-contraction, 
too), we finally obtain 

||(Zn+l ~ llc“(R‘^,R) — 11^ ~ 9"-l) llc“(R‘^,R) — ll?" ~ Qn-1 llcO(R'^,R) ' 
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□ 


Example 1 (Bermudan put option with equidistant exercise times 
in t ■ No on the weighted arithmetic average of a basket in a discrete 
Markov model with a discount factor c = e~'~* for r > 0). Let /3i,Pd S 
[0,1] be a convex combination and assume that A is such that 


m 



( 2 ) 


then the functions 


d 



and h •= K (where K > 0) satisfy the equations Ah = h and Ag = g, re¬ 
spectively. Moreover, by definition g < h. Then we know that the (perpetual) 
Bermudan option pricing algorithm that iteratively applies V to the payoff func¬ 
tion gV 0 on the log-price space, will increase monotonely and will have a limit 
which is the smallest nonnegative fixed point of V. Moreover, the convergence 
is linear and the contraction rate can be bounded by c. 

The condition ^ can be achieved by a change of the time scale (which 
ultimately leads to different cubature points for the distribution of the asset price) 
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